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ABSTRACT 

Wc present weak gravitational lensing analysis of 22 high-redshift (z > 1) clusters based on Hubble 
Space Telescope images. Most clusters in our sample provide significant lensing signals and are well 
detected in their reconstructed two-dimensional mass maps. Combining the current results and our 
previous weak-lensing studies of five other high-2: clusters, we compare gravitational lensing masses 
of these clusters with other observables. We revisit the question whether the presence of the most 
massive clusters in our sample is in tension with the current ACDM structure formation paradigm. 
We find that the lensing masses are tightly correlated with the gas temperatures and establish, for 
the first time, the lensing mass-temperature relation at 2; > 1. For the power law slope of the M — Tx 
relation (M oc T"), we obtain a = 1.54 ± 0.23. This is consistent with the theoretical self-similar 
prediction o; = 3/2 and with the results previously reported in the literature for much lower redshift 
samples. However, our normalization is lower than the previous results by 20-30%, indicating that 
the normalization in the M — Tx relation might evolve. After correcting for Eddington bias and 
updating the discovery area with a more conservative choice, we find that the existence of the most 
massive clusters in our sample still provides a tension with the current ACDM model. The combined 
probability of finding the four most massive clusters in this sample after the marginalization over 
cosmological parameters is less than 1%. 

Subject headings: gravitational lensing — dark matter — cosmology: observations — X-rays: galaxies: 
clusters — galaxies: clusters — galaxies: high-redshift 
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Gravitationally bound systems have been the main 
subjects of studies for astronomers throughout history 
because they are observationally identified as discrete 
entities in the universe. Galaxy clusters, the largest 
among these, are believed to be also the last systems 
detached from the Hubble expansion. Sttidying galaxy 
clusters provides unique opportunities to test our struc- 
ture formation paradigm, understand gas physics, infer 
the properties of dark matter, investigate cluster galaxy 
evolution, and constrain cosmological parameters. 

The formalism for the last of these, i.e., the use of clus- 
ters for cosmology, was pioneered by Press & Schechter 
(1974), who realized that under the assumption of Gaus- 
sian primordial density fluctuation the fraction of mass 
in halos more massive than the threshold M is related 
to the fraction of the volume in which the smoothed ini- 
tial density field is above some threshold density contrast 
6c- This simple, but ingenious, approach has been shown 
to work well in comparison with numerical simulations. 
Since this early work, many authors have extended the 
formalism to provide better agreement with the recent 
state-of-the-art numerical simulations. 

The number of clusters of a given mass, or the mass 
function, in the local universe constrains the combined 
properties of the matter density and its fiuctuation 
^Mcrl^ (e.g., Bahcall & Cen 1993; PierpaoH et al. 2001; 
Reiprich & Boehringer 2002, and references therein) 
whereas the evolution of the mass function breaks this 
degeneracy (e.g., Bahcall & Fan 1998). Therefore, there 
has been continuous effort to extend our knowledge of 
the cluster mass function to higher and higher redshift 
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regime. During the last decade many on-going surveys 
have increased the number of known of z > 1 clusters by 
many factors (e.g., Eisenhardt et al. 2008; Boehringer et 
al. 2006; Gladders & Yee 2005). Therefore, it is now pos- 
sible to study the evolution of the mass function across 
a large range of redshifts. Such a study can not only 
provide an important check to the results so far obtained 
from lower redshift samples, but also enable a critical 
test on the primordial non-gaussianity (e.g., Jimenez & 
Verde 2009; Sartoris et al. 2010). 

Without question, one of the most crucial issues is the 
accurate derivation of masses for high- redshift clusters. 
Indirect approaches, i.e., the estimation from velocity 
dispersion. X-ray properties, or Sunyaev-Zeldovich ob- 
servations, require assumptions on the dynamical state of 
the clusters and/or the mass versus mass proxy calibra- 
tion. By contrast, gravitational lensing is a unique tool to 
obtain the cluster mass without relying on any dynamical 
assumption. However, because this method requires ex- 
pensive, high-resolution observations from space, it is not 
practical, if not impossible, to perform a lensing analysis 
for complete samples. Instead, it is much more feasible 
to apply the technique to a subset of the sample, and to 
use the results to calibrate the relation between masses 
and their proxies. 

In this paper, we present weak-lensing analysis of 22 
z > 1 clusters based on Hubble Space Telescope images 
(HST). Our detailed study of these clusters via their 
weak gravitational lensing signal will provide the afore- 
mentioned, invaluable calibration between weak-lensing 
masses and other observables. In addition, the most 
massive among these allow us to investigate the behav- 
ior of the high end of the mass function (Hoyle et al. 
2011; Holz & Perlmutter 2010). The 22 clusters were 
discovered in different surveys, and do not represent a 
complete sample, implying that any significant excess of 
massive clusters (beyond the sample fluctuation) would 
in fact underestimate the difference between observation 
and theory. In order to complement our sample for the 
investigation of the mass and mass observable relation 
and the implication for the cosmology, we include five 
other high-rcdshift clusters that were previously studied 
via lensing. 

The structure of the current paper is as follows. In 
S|2] we describe our ACS data. The method and the 
weak lensing mass reconsruction for the 22 clusters are 
presented in S|3j S|4] compares the lensing measurements 
with other cluster properties for the combined sample of 
27 clusters and presents the mass-temperature relation at 
z > 1. We investigate the abundance of the most mas- 
sive clusters in our study relative to the theoretical ex- 
pectations within the current paradigm of the structure 
formation in f|5] before the conclusion in S|6j Discussions 
on the details of our corrections to shape measurements 
including point spread function (PSF) and charge trans- 
fer inefficiency (CTI) are deferred to Appendices A and 
B, respectively. 

We assume (f^Mj^^A,^) — (0.3,0.7,0.7) for cosmology 
unless explicitly stated otherwise. All the quoted uncer- 
tainties are at the 1-cr (~ 68%) level. 

2. OBSERVATIONS 

The 22 z > 1 clusters were observed with ACS/WFC 
using the F775W and F850LP filters (hereafter referred 



to 1775 and Z850, respectively) as part of the HST Clus- 
ter Supernova Survey (GO-10496, PI: Perlmutter) during 
the period of 2005 July-2006 November (see Dawson et 
al. 2009). Three clusters among these were also observed 
as part of Guaranteed Time Observation 9290 and 9919 
(PI: Ford) during the periods of 2002 May- June and 2004 
February-March. Table 1 summarizes the exposure times 
of the 22 clusters in ■8775 and zsso- 

2.1. ACS Data Reduction 

Our reduction started with the FLT images, which are 
the products of the standard STScI CALACS pipeline 
(Hack et al. 2003). We used the "apsis" pipeline 
(Blakeslee et al. 2003) to determine shift and rotation, 
correct geometric distortion, remove cosmic rays, and 
create final mosaic images. The determination of the 
shift and rotation between different pointings is the most 
critical step among these for a weak-lensing analysis be- 
cause any potential misalignment induces a coherent dis- 
tortion of object shapes, mimicking gravitational lensing 
signals. The "apsis" pipeline uses high S/N astronomical 
objects iteratively to measure and refine the alignment 
through the "match" program developed by Richmond 
(2002). We find that the average uncertainties for shift 
and rotation are ~ 0.02 pixels and ~ 0.2", respectively, 
and thus the systematics, if any, due to the alignment er- 
rors is negligible. The LanczosS (windowed sine function) 
kernel was used for drizzling (Fruchter & Hook 2002) 
with the native 0.05" pixel scale of ACS/WFC. 

2.2. Object Detection and Source Selection 

A detection image was created for each cluster by 
weight-averaging the two passband images via apsis. 
Then, SExtractor (Bertin & Arnouts 1996) was run in 
dual image mode so that objects were identified from 
the detection image while photometry was measured on 
an individual filter. This provides common isophotal ar- 
eas to both filters, which enables a robust object color 
measurement. We filtered the detection image with a 
Gaussian kernel that approximately matches the PSF of 
the instrument before looking for 5 or more contiguous 
pixels above 1.5 times the sky rms. Although our crite- 
ria were experimentally determined to minimize the re- 
turn of false detections by the software, inevitably man- 
ual identification of spurious detections (e.g., saturated 
stars, diffraction spikes, uncleaned cosmic rays near field 
boundaries, H II regions inside nearby galaxies, etc.) is 
always required. We divided the 22 clusters into several 
groups and a few authors were assigned to each group to 
carry out this manual procedure. A possible concern is 
whether or not different groups might have non-negligible 
biases in this false object identification, which may prop- 
agate into cluster mass determination. To examine the 
possibility of bias in the manual procedure, several clus- 
ters were randomly selected by one author, and com- 
pared the two false object catalogs created by this author 
and others. We find that about ~ 60% of false objects 
are shared by the two catalogs, and the disagreement is 
mostly attributed to very faint objects, whose fluxes are 
just above the sky background. Because these extremely 
faint objects were discarded anyway by our magnitude 
and shape measurement error criteria, the resulting dif- 
ference in mass estimation and two-dimensional mass re- 
construction is far below statistical errors. 
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We define source galaxies as objects whose 1775 — zgso 
colors are bluer than the color of the red-sequence in each 
cluster while their total ^sso magnitudes (MAG_AUTO) are 
fainter than zsso = 24. We also require the source galax- 
ies to have ellipticity measurement error less than 0.25. 
The total exposure times for the 22 clusters vary sub- 
stantially and thus so does the number density of source 
galaxies as shown in Table 1. 

2.3. Shape Measurement 

Our shape measurement method is detailed in Jee et 
al. 2009. Briefly, we fit a PSF-convolved elliptical Gaus- 
sian function to source galaxies to determine their semi- 
major and -minor axes. Convolution with an elliptical 
Gaussian is required to measure galaxies' ellipticity be- 
fore photons reach the instrument. Because an elliptical 
Gaussian profile is not an unbiased representation of true 
galaxies, the measurement inevitably introduces a bias. 
On average, the bias is toward underestimation of ellip- 
ticity because a Gaussian profile is steeper than that of 
real objects; more extended objects are subject to larger 
underestimates of ellipticity. In future analyses, it may 
be worth experimenting with more generalized functional 
forms and examining the effect based on simulated im- 
ages. Nevertheless, we stress that the amount of bias 
induced by the current shape measurement is small and 
to first order was corrected here using the simulation re- 
sults of Jee et al. (2007b), which shows a ^ 6% bias for 
7^ 0.5. 

We model the PSF using our library constructed from 
stellar field observations (Jee et al. 2007b). In order to 
find the matching PSF template, we use ^ 10 high S/N 
stars present in each cluster field. This matching is done 
for each visit, and the final PSF model is obtained after 
the model for each visit is shifted, rotated, and stacked. 
The mean residual ellipticity is less than ~ 0.01, and 
this level of accuracy is sufficient for the current cluster 
lensing analysis. The detailed results for each cluster can 
be found in Appendix \K\ 

Together with the PSF, CTI may be a potential source 
of systematics in weak lensing. For bright stars, the trail- 
ing from deferred charges is visually apparent. We quan- 
tify the effect using cosmic rays and warm pixels, and 
find that the effect on the ellipticity of the galaxies that 
we use to extract the signal is in fact much less than 
the statistical noise set by the finite number of galaxies. 
Thus, we confirm our previous argument (Jee et al. 2009) 
that the ACS CTI effect is negligible for cluster lensing 
analysis although it is a critical source of systematics in 
cosmic shear studies. We present the details in Appendix 

m 

3. WEAK-LENSING ANALYSIS 

3.1. Mass Estimation and Two-dimensional Mass 
Reconstruction 

In the weak-lensing regime, the characteristic length of 
the lensing signal variation is larger than the object size. 
Thus, the resulting shape distortion can be linearized as 
follows: 



A(x) = 4 



a25f(x) 

dxidxi 



1 - K - 71 



-72 

K + 7i 



where A(x) is the transformation matrix x' = Ax, which 
relates a position x in source plane to a position x' in 
image plane, and ^ is the two-dimensional lensing po- 
tential. The convergence k is the surface mass density in 
units of critical surface mass density 
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where D{zs), D{zi), and D{zi,Zs) are the angular diam- 
eter distance from the observer to the source, from the 
observer to the lens, and from the lens to the source, re- 
spectively. The convergence k, and the shears 71(2) are 
related to the lensing potential ^ via 

^ = ^(^11 + ■^22), 71 = ^(^11 ~ ■^22), and 72 = i/'i2 = V'^Si) 

where the subscripts on V'i(j) denote partial differentia- 
tion with respect to Xi(^jy 

A useful quantify to estimate masses is the reduced 
tangential shear defined by 
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where <j) is the position angle of the object with respect to 
the cluster center, and g is the reduced shear 7/(1 — n) 
(valid only in the weak-lensing regime). This reduced 
tangential shear profile informs us of how the mass den- 
sity of the lens changes as a function of radius. Of 
course, without any lensing signal, the resulting shear 
profile should vanish with ffuctuations consistent with 
shot noise. 

Hoekstra (2001; 2003) demonstrates that a cosmic 
shear (lensing by large scale structure) is an important 
limiting factor in the accuracy of cluster masses derived 
from the cluster tangential shears. Following the for- 
malism of Hoekstra (2003), we estimate that the cosmic 
shear contribution on average is 7 ~ 0.01 for our sample. 
Although this is small compared to the shot noise set by 
the finite number of source galaxies, we include the effect 
in our final mass uncertainty. 

A potential ambiguity is the choice of the cluster cen- 
ter when there is disagreement between centroids deter- 
mined from the distribution of cluster galaxies, lensing 
mass peaks, and X-ray emission. If we adopt a location 
that maximizes the amplitude of the tangential shear or 
a peak in the convergence map as the cluster center, the 
derived mass will be always biased high. Fortunately, in 
previous high S/N weak-lensing studies, where the num- 
ber density of background galaxies is high (> 120 per sq. 
arcmin), the weak-lensing mass peaks are in good spa- 
tial agreement with cluster galaxies. However, in many 
cases the X-ray peaks are conspicuously offset, indica- 
tive of merging activity or the collisional properties of 
the intracluster medium. In this paper we adopt the 
centroid defined by the cluster galaxies as the center for 
the construction of the tangential shear profile. For the 
clusters in our sample where the statistical significance of 
the lensing signal is high, these centroids also agree well 
with the mass peaks. In several cases, where the lens- 
ing signal is weak (the cluster is not massive and/or the 
image is not sufficiently deep), we occasionally observe 
rather large (> 20") offsets between mass and galaxies. 
The choice of centroid will lead to an underestimate of 
the cluster mass if the offsets refiect real features in the 
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system. We do not use the inner most data points in our 
mass estimation because the tangential shears at very 
small radii are highly sensitive to the centroid choice. 
Furthermore, both observationally and theoretically, it 
is not clear how the cluster mass profile behaves near the 
center. 

There are two popular ways in deriving the lensing 
mass from tangential shear: aperture-mass densitometry 
(Fahlman et al. 1994) and parametrized model fitting. 
Aperture-mass densitometry is useful when one attempts 
to estimate a total projected mass within some aperture 
radius without requiring an assumption on the behav- 
ior of the cluster mass profile. However, this approach is 
not practical for the current data set, which in most cases 
provides only 3' x 3' areas smaller than the virial radii 
of the clusters. Therefore, we use the second method to 
derive the cluster masses. This method determines the 
parameters of analytic halo models by comparing the ob- 
served tangential shears with the expected values. Of 
course, a scatter is introduced because the real cluster 
profile is different from the model. Nevertheless, many 
studies show that this parametric approach gives consis- 
tent results with the values obtained by the former (e.g., 
Jee et al. 2005a; Hoekstra 2007). 

We employ two kinds of halo models: singular isother- 
mal sphere (SIS) and Navarro-Frenk- White (NFW; 
Navarro et al. 1997) profiles. Although the SIS profile is 
considered inappropriate at large radii (many lines of ev- 
idence suggest the mass density drops faster than l/r^), 
the resulting parameter can be conveniently translated 
to the Einstein radius Oe ov the velocity dispersion of 
the system. We use this predicted velocity dispersion to 
compare with the dynamical value. When we fit NFW 
profiles, it is assumed that the cluster virial mass M2Q0 
( the total mass at the radius inside of which the mean 
density is 200 times the critical density of the universe 
at the cluster redshift) is tied to the concentration c via 
the following relation of Duffy et al. (2008): 



5.71 



2 X 1012 /i-i Mr, 



(1 + ^) 



-0.47 



(5) 



Therefore, effectively, our NFW model is described by a 
single parameter. The relation between projected mass 
density and observed shear is simple for SIS with an Ein- 
stein radius rg: k = 0.5r£:/r and g = k/{1 — k) (in the 
weak-lensing regime). For NFW, the relation is rather 
complicated, and we refer readers to Bartelmann (1996). 
We present tangential shears and fitting results for indi- 
vidual clusters in Figure JT][2] 

In principle, the two-dimensional projected mass dis- 
tribution K can be obtained by convolving the shear 7 as 
follows (Kaiser & Squires 1993): 



«(x) = i J D*{y 



c')7(x')d'x. 



(6) 



where D{x.) = — l/(xi — 1x2)^ is the convolution ker- 
nel. However, the direct application of Equation [6] gives 
rises to some practical problems. First, the ellipticity of 
individual galaxies should be smoothed, and we do not 
know the optimal smoothing scale beforehand. Second, 
the smoothed galaxy ellipticity gives only an estimate 
for the reduced shear g = 7/(1 — k). Third, the rela- 
tion 5 = 7/(1 — k) is only valid in the [7/(1 — k)\ < 1 
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regime. Fourth, it is not obvious how to weight shape 
measurement noise properly via Equation [6] 

Certainly, a more robust two-dimensional mass recon- 
struction algorithm is required to account for the afore- 
mentioned subtleties and minimize numerous artifacts. 
Many suggestions are present in the literature (Bridle et 
al. 1998; Seitz et al. 1998), and the consensus is that a 
reliable mass reconstruction should be achieved through 
iterations. In addition, it is desirable to use the ellipticity 
of individual galaxies rather than the smoothed values 
so as to reveal small-scale, but significant features. In 
the current paper, we reconstruct two-dimensional mass 
maps using the entropy-regularized, maximum likelihood 
code of Jee et al. (2007a). This method of mass re- 
construction accounts for the aforementioned subtleties 
and is effective in revealing low-contrast str uctu res. We 

, , Results 

on individual clusters are discussed in 

3.2. Redshift Estimation of Source Galaxies 

Because Sc scales as oc D{zs)/[D{zi)D{zi, Zs)] (eqn. 
[2|, the mass density unit rises sharply as the redshift 
of the lens approaches that of the source. Thus, in the 
weak-lensing analysis of high-redshift clusters, small er- 
rors in our estimate of the source redshift distribution 
result in large errors in the translation of the signal into 
the physical unit. For example, a 10% systematic uncer- 
tainty in the effective source redshift leads to 1 — 2% error 
in mass for a z = 0.2 object (e.g., ABELL 1689). The 
same amount of uncertainty would give ~ 30 % error for 
the mass of a z = 1.4 lens (e.g., XMM2235-2557Q 

We estimate the redshift distribution of the source 
population using the publicly available Ultra Deep Field 
(UDF) photometric redshift catalog (Coe et al. 2006). 
The ultra deep images in six filters from F435W to 
F160W provide unparalleled high-quality photometric 
redshift information well beyond the limiting magnitude 
of the cluster images that we analyze here. The catalog 
has been extensively used in our previous studies (e.g., 
Jee et al. 2009), and thus we only briefly summarize the 
procedure. 

First, we bin our source galaxies in 0.5 mag intervals. 
Next, we randomly draw the galaxies that match the se- 
lection criteria of each bin from the UDF photo-z catalog. 
Finally, we combine these simulated catalogs for all bins 
to create the redshift catalog for the entire source popula- 
tion. Because the UDF is much deeper than the cluster 
images and also because there is sample variance, it is 
important to take into account the difference in number 
density in this step. Consequently, we are utilizing the 
magnitude-color-redshift relation measured in the UDF 
data, rather than simplistically imposing the UDF red- 
shift distribution on our images, which would cause much 
greater systematics. 

As noted in Jee et al. (2009), the total error in Zg// 
due to he sample variance, the resampling error, and 
the difference among the photo-z estimation codes is 
Szeff ~ 0.06. This causes ^ 11% (^ 3 %) uncertainty in 
mass for a z ^ 1.4 (z ~ 0.9) cluster. We include these er- 
rors in our final uncertainty. The assumption that there 

The exact value also depends on the depth of the image, 
which determines how many distant galaxies are resolved. Here we 
assume ~ 4 orbit integration with HST/ACS. 
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exists a single source plane while in reality each source 
galaxy is at different redshift also creates a bias in the in- 
terpretation of the lensing signal. We correct this bias to 
first order using the equation derived by Seitz & Schnei- 
der (1997). 

3.3. Notes on Individual Clusters 

We here comment on specific features in the measure- 
ment of each cluster. All of the key numbers are sum- 
marized in Table 2. 

3.3.1. XMMXCS J2215-1738 (z = 1.46; 

XMMXCS2215-1738 was the highest redshift cluster 
spectroscopically confirmed (Stanford et al. 2006) un- 
til the recent discovery of CIG J0218. 3-0510 at z = 
1.62 (Papovich et al. 2010). Assuming no signifi- 
cant point source contamination in the archival XMM- 
Newton data, Stanford et al. (2006) obtained a temper- 
ature of Tx = 7.4^^ g keV with a bolometric luminosity 
of Lx = 4.4^°;^ X lO^^ergs s'^. Hilton et al. (2010) 
reported significantly lower values Tx — 4.1^Q g keV and 

Lx = 2.92^0 35 X lO'^'^ergs based on a joint analysis of 
Chandra and XMM-Newton data, which enables them 
to remove the point source contamination from the dif- 
fuse component. Assuming an isothermal /3 model with 
(3 = 0.63 and = 52 kpc, the temperature Tx = 4.1^Qg 

keV translates into M200 = S.Otg g x 10" M©. 

We initially considered a weak-lensing study of 
XMMXCS2215-1738 to be highly chaUenging because of 
the low X-ray mass, high redshift, and relatively shallow 
depth of the F775W image, where we measure galaxy 
shapes. However, the cluster is clearly detected in our 
two-dimensional mass reconstruction, which shows the 
mass centroid in good spatial agreement with the cluster 
members and the X-ray centroid (Figure pi) . The weak- 
lensing mass of the cluster is M200 = 4.31^^ ^ 10^^ -^O; 
a factor of two higher than the X-ray estimate (however, 
both results are marginally consistent). 

The predicted velocity dispersion 9421; ijje km is 
consistent with the spectroscopically measured value 
710 ± 110 km derived from 31 members within i?200 
(Hilton et a. 2010); interestingly, if all 44 members 
were used, the velocity dispersion increases to 890 ± 
110 km s~^ giving a better agreement with the lensing 
result. 

3.3.2. XMMU J2205-0159 (z = 1.12; 

XMMU2205-0159 was discovered in the XDCP 
(Boehringer et al. 2006) survey in the archival image 
of QSO B2202-0209 at z=1.77. The mass centroid is 
offset north of the BCG (a, ^) -(22:05:50.7,-01:59:30) 
by ^ 17" (Figure [4]). We derive a weak-lensing mass of 

^^200 — S.Oii^o ^ -^0> making the cluster among the 
least massive clusters in our sample. Neither a dynami- 
cal nor X-ray study of the cluster has been published to 
date. 

3.3.3. XMMU ,71229+0151 (z^ 0.98) 

XMMU J1229+0151 was serendipitously detected in 
the field targeted for the well-known quasar 3C 273. 
The cluster is hot and X-ray luminous with Tx — 



6Atal keV and = 3 x 10''^ erg s"i (Santos et 
al. 2009). Our lensing analysis predicts a velocity dis- 
persion of 867lgg km s^^ and gives a virial mass of 
-^200 = 5.3l^'2 X lO^'' Mq. This lensing mass agrees 
nicely with the X-ray mass M200 = 5.7tJ;8 x lO" Mq 
when an isothermal (3 model with /3 = 0.7 and Tc = 250 
kpc is assumed; no measurement of the X-ray surface 
brightness has been reported yet. However, the velocity 
dispersion 683±62 km s""'^ derived from 27 spectroscopic 
members (Santos et al. 2009) is significantly lower than 
the lensing prediction 867^69 km s^^. Our convergence 
map shows a strong peak ^ 10" northeast of the X-ray 
peak and the cluster galaxies (Figure [5]) . 

3.3.4. WARPS J1415+3612 (z = 1.03^ 

The cluster WARPS J1415-f-3612 was discovered (Perl- 
man et al. 2002) in the first phase of the Wide Angle 
ROSAT Pointed Survey (WARPS). The ICM tempera- 
ture of the cluster has been measured from the Chandra 
(5.59 ± 0.84 keV; Allen et al. 2008) and XMM-Newton 
(5.7j:J:7 keV; Maughan et al. 2006) data analysis. The 
cluster has been considered relaxed because of its sym- 
metric X-ray emission (Maughan et al. 2006; Allen et 
al. 2008). Our weak-lensing mass map of the cluster also 
does not show any significant substructure, adding an ad- 
ditional support to the hypothesis. Both virial masses de- 
rived from X-ray and lensing (M200 = 4.6to'8 ^ ^^^^ -^o 
and 4.71^4 ^ ^^^"^ Mq, respectively ) are in excellent 
agreement. 

The ACS image of the cluster shows a strongly lensed 
arc at z ~ 3.898. The source is a Ly-a emitter, and is 
located ~ 7" away from the BCG (Huang et al. 2009). 
When we adopt this distance as the Einstein radius at 
z = 3.898, the value is consistent with the weak-lensing 
prediction Bsiz = 3.898) = 8" ± 2" derived from the 
NFW fitting result; note also that the location of the 
weak-lensing mass peak is in good agreement with that 
of the BCG (Figure |6|. 

3.3.5. ISCS JI432+3332 (z = 1.11; 

The ISCS J1432-K3332 cluster was reported by Elston 
et al. (2006) as one of the first spectroscopically con- 
firmed 2 > 1 clusters detected through the joint use of 
the FLAMINGOS Extragalactic Survey (FLAMEX; El- 
ston et al. 2006), the NOAO Deep Wide-Field Survey 
(NDWFS; Brown et al. 2003), and the Spitzer IRAC 
Shallow Survey (ISCS; Eisenhardt et al. 2004) data sets. 
They estimated from 8 spectroscopic members that the 
velocity dispersion of the cluster is 920 ± 230km s~^. 
Eisenhardt et al. (2008) reported a refined measure- 
ment of CFy = 734 km s~^ using 23 redshifts; the au- 
thors did not quote an uncertainty and we estimate 
5a ^ ^ 115 km sT^ . The virial mass is estimated to be 
M200 = 4.9li[;2 X 10^^ Mq- The mass contours are weh 
traced by the cluster galaxies (Figure [t]). 

3.3.6. ISCS JI429+3437 (z = 1.26) 

ISCS J1429-I-3437 was identified photometrically with 
the combined use of the ISCS and NDWFS data sets 
(Eisenhardt et al. 2008). From the 9 spectroscopic mem- 
bers of Eisenhardt et al. (2008), we calculate the velocity 
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dispersion to be 767±295 km s"\ which agrees well with from our lensing analysis, which give 833ii5o km 



the predicted velocity dispersion 732^7° km s ^ . For the 



virial mass, we obtain Af2oo = 5.4^^ g x 10^* Mq. 

We detect two mass clumps within the ACS field of the 
cluster observation (Figure |8|. The western mass peak is 
much stronger and spatially correlated with the cluster 
galaxy candidates. The much weaker eastern clump does 
not have any compact galaxy distribution. Nevertheless, 
the global east-west elongation of the mass distribution 
seems to follow the early-type galaxies in the field whose 
colors are consistent with that of the red-sequence at the 
redshift of the cluster. 

3.3.7. ISCS JI434+3427 (z = 1.24) 

The ISCS J1434-h3427 cluster, discovered in the ISCS 
and NDWFS survey, is reported to possess a pronounced 
filamentary structure (Brodwin et al. 2006). Although 
our mass reconstruction reveals a significant mass peak 

10 " east of the compact galaxy distribution, the current 
1775 image is not deep enough {'^ 2685 s) to study the 
substructure of this high-redshift cluster in detail. Be- 
cause the cluster is located near the edge of the ACS 
field, we need to assume an azimuthal symmetry at large 
radius to obtain a mass estimate. This may substantially 
bias our measurement if the cluster departs significantly 
from the assumed axisymmetry. Nevertheless, we note 
that the predicted velocity dispersion 770^^33 km s~^ by 
our lensing analysis is consistent with the dynamical ve- 
locity dispersion 863 ± 170 km s~^ that are derived from 

11 cluster members (Meyers et al. 2011). 

3.3.8. ISCS JI432+3436 (z = 1.35) 

The cluster was discovered in the ISCS and NDWFS 
survey, and about 8 members have been spectroscopi- 
cally confirmed (Eisenhardt et al. 2008). The cluster 
detection is very strong in our mass reconstruction de- 
spite both the high-redshift and the shallowness (2235 
s) of the «755 image. We estimate the virial mass to be 
M200 = X 10" Mq. The velocity dispersion de- 

rived from 8 members is 807 ± 340 km s~^, which agrees 
well with the lensing prediction 9121^'^q2 km s^^. 

3.3.9. ISCS JI434+3519 (z = 1.37; 

Eisenhardt et al. (2008) report that five members have 
been spectroscopically confir med . Our two-dimensional 
mass reconstruction (Figure 11) reveals a weak detec- 
tion of convergence peak near the cluster center defined 
by Eisenhardt et al. (2008) using photometric redshifts. 
The 15" offset may be attributed to statistical noise be- 
cause of the high redshift, the low mass, and the insuf- 
ficient depth (1920 s) of the 1775 image. The virial mass 
is determined to be M200 = 2.8ti;4 x lO" ^Q- 

3.3.10. ISCS JI438+3414 (z = 1.41; 

The discovery of ISCS J1438-h3414 is reported by Stan- 
ford et al. (2005), who confirmed five spectroscopic mem- 
bers. Brodwin et al. (2011) quote a dynamical velocity 
dispersion of 757^203 from a total of 11 cluster 

members, which can be converted to A/200 — 2.3l2'^ ^ 
10^^ Mq using the mass-dispersion relation from Evrard 
et al. (2008). These values are consistent with those 



and a virial mass of A/200 = 3. 1^1:4 x lO^"* Mq. An- 
dreon et al. (2011) measured an X-ray temperature 
Tx = 4.9t?:6 keV from the relatively deep (~ 150 ks) 
Chandra data of the cluster. This temperature is con- 
sistent with the measurement (3.3ii[ g keV) of Brodwin 
et al. (2011) obtained from the same data. Assuming 
P = 0.7 and rc - 100 kpc, we convert Tx = 4.9t5 j keV 

to A/^Q = 3.2^^ 4 Mq, which also agrees well with our 
lensing mass. The centroid of the diffuse X-ray emission 
is in good spatial agreement with that of the weak-lensing 
mass. 

3.3.11. RCS 0220-0333 (z = 1.03^ 

The clusters with a prefix RCS hereafter were mostly 
discovered in the Red-sequence Cluster Survey-I (Glad- 
ders & Yee 2005); the exception is RCS 1511-^0903, 
which was discovered in the Red-sequence Cluster 
Survey-II (Gilbank et al. 2011). RCS 0220-0333 is 
an optically rich cluster at z = 1.03. In the ACS 
pseudo-color composite image, the early- type galaxies 
with 1775 — ^850 ~ 0.9 appear to form a North-South fil- 
amentary structure. The cluster is strongly detected in 
lensing. The mass centroid lies close to the cluster mem- 
ber with pronounced strong lensing features. However, 
the mass distribution does not indicate the North-South 
elongation seen in the cluster galaxies. We obtain an 
virial mass of M200 = 4.8j;^;3 x lO"^'' Mq with a predicted 
velocity dispersion of 881]'^74 km s^^. 

3.3.12. RCS 0221-0321 (z = 1.02) 

This optically rich galaxy cluster RCS 0221-0321 is 
clearly visible in our mass reconstruction. The mass 
map show a significant mass clump aligned with the 
optical center (e.g., see the isodensity contours for red 
galaxies from Andreon et al. 2008). However, this 
cluster is found to be one of the least massive clus- 
ters in our sample. We estimate the virial mass to be 
A/200 = l-8to.? X 10^^ ^0- The predicted velocity dis- 
persion 699^94 km s~^ agrees well with the dynamical 
measurement 710±150 km s~^ based on 21 spectroscopic 
members (Andreon et al. 2008). 

3.3.13. RCS 0337-2844 (z = 1.10) 

RCS 0337-2844 is clearly visible in our two-dimensional 
mass map, which shows a good mass-galaxy correlation. 
The exposure time of the 1755 image of this cluster is 
the shortest among the clusters studied here, which in- 
dicates that this strong weak-lensing signal is due to a 
significant mass. We estimate A/200 = 4.9l^'7 ^ ^^^^ ^Q- 
No spectroscopic data have been published to date. 



3.3.14. RCS 0439-2904 (z = 0.95) 

Barrientos et al. (2004) reported that RCS 0439-2904 
is an optically rich cluster (~ 9 a in the Ks image) at 
a high redshift. Cain et al. (2008) measured the X-ray 
temperature of the cluster from the 200ks Chandra to be 
Tx = 1.5^Q 4 keV, which is significantly lower than what 
the optical richness suggested. The single isothermal /3 
model gives A/200 = ^-Qtl " x 10" Mq. On the other 
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hand, the mass-richness relation implies a much higher 
mass by at least an order of magnitude. Because their 
velocity histogram show that there are multiple compo- 
nents, Cain et al. (2008) interpreted this large discrep- 
ancy as indicating the presence of two clusters along the 
line-of-sight direction. They also claim that the unex- 
pectedly high gas fraction obtained when a single com- 
ponent is assumed supports the two-component interpre- 
tation. 

The line-of-sight hypothesis makes our lensing study 
of the cluster interesting. If there is indeed a line-of- 
sight superposition of two clusters, the mass estimate 
from lensing should give the sum of the two clusters 
as opposed to the single component X-ray model. Our 
analysis shows that the cluster virial mass is M20Q — 
4.3t X 10^"* A/0 , nearly a factor of 10 higher than the X- 
ray prediction as was also indicated by the mass-richness 
relation. Interestingly, our predicted velocity dispersion 



831 



-1-68 
-74 



km s is consistent with the dynamical mea- 
surement 1080 ± 320 km ■sT^ of Cain et al. (2008). The 
lensing mass is consistent with the mass-richness relation 
if the cluster is assumed to consist of two clusters; the 
new lensing mass shifts the Model II data point of Figure 
3 of Cain et al. (2008) upward in such a way that the 
new point is well bracketed by the 1 — a scatter of the 
mass-richness relation of Blindert (2006). 

3.3.15. RCS 2156-0448 (z = 1.07; 

The cluster RCS 2156-0448 was reported as an opti- 
cally rich cluster with a strong lensing arc candidate by 
Gladders et al. (2003). If confirmed, the presence of 
this arc suggests that the cluster is massive as indicated 
by its optical richness. However, Gladders et al. (2003) 
commented that this candidate should be considered the 
least likely among their secondary sample. 

Our weak-lensing mass reconstruction reveals only a 
weak convergence peak near the arc candidate. Instead, 
the mass map shows a more significant mass peak ^ 70" 
south of the assumed cluster center. Our visual inspec- 
tion of the ACS image shows that on this location there 
seems to be also early-type galaxies whose colors are con- 
sistent with that of the cluster red-sequence. 

Nevertheless, one must be reminded that the combined 
(i775 + -zsso) ACS image of the cluster is the shallow- 
est among our 22 cluster sample, which gives the lowest 
number density of background galaxies (~ 50 arcmin^^). 
Hence, the mass-galaxy comparison should await a future 
analysis with deeper images. Our mass presented here 
is estimated by placing the center on the strong-lensing 
candidate of Gladders et al. (2003). 

3.3.16. RCS 1511+0903 (z = 0.97^ 

The ACS image of RCS 1511-K0903 shows a compact 
distribution of early-type galaxies whose colors are con- 
sistent with the cluster redshift z = 0.97. The cluster 
is clearly visible in our mass reconstruction, which re- 
veals a strong mass peak at the location of the cluster 
galaxies. Based on 9 redshifts, Meyers et al. (2011) 
estimate the velocity dispersion of the cluster to be 
717 ± 208 km s~^, which is consistent with our lensing 
prediction 699;tio9 "^^^ virial mass of the cluster 

i.9tJ:^ X 10^4 Mq. 



from our lensing analysis is M; 



200 



3.3.17. RCS 2345-3632 (z = 1.04; 

RCS 2345-3632 is an optically rich cluster at z=1.04 
with 23 spectroscopically confirmed cluster galaxies. The 
dynamical velocity dispersion 670 ± 190 km sT^ is con- 
sistent with the lensing prediction 684j^yg km s^^. The 
two-dimensional mass map agrees well with the cluster 
galaxy distribution (Figure 19). The relatively deep ACS 
image allows us to utilize 103 background galaxies per 
sq. arcmin, and thus the resulting mass reconstruction 
should be considered most reliable among the 8 RCS clus- 
ters presented in this paper. We estimate that the virial 
mass of the cluster is M200 = x 10^"* Mq. 

3.3.18. RCS 2319+0038 (z = 0.9i; 

RCS 2319-1-0038 is a strong-lensing cluster showing at 
least two spectacular arcs (Gladders et al. 2003). The 
recent ACS images of the cluster reveal more strong- 
lensing features (at least 7 tangential and 2 radial arc 
candidates). Gilbank et al. (2008) discovered that the 
cluster is in fact part of a supercluster containing two 
other clusters RCS 2319+0030 and RCS 2320+0033. All 
three clusters are well detected in the Chandra X-ray ob- 
servations with similar gas temperatures ~ 6 keV (Hicks 
et al. 2008). RCS 2319+0038 is separated from RCS 
2320+0033 and RCS 2319+0030 by - 5' and - 7', re- 
spectively, and therefore our lensing analysis based on 
the ACS image covering the central 3' x 3' region of RCS 
2319+0038 is not likely to be affected by these two other 
structures. Gilbank et al. (2008) reported a velocity dis- 
persion of 990 ± 240 km s~^, which agrees well with the 
current lensing prediction 898^^} km Our conver- 
gence map reveals a single strong mass peak at the center 
of the strong-lensing system without any significant sub- 
structures. Both this mass map and the tangential shear 
profile suggest that the cluster is relaxed. The virial 
mass derived from the X-ray measurements of Hicks et al. 
(2008) {Tx = 6.2^°;^ keV, = 100 kpc, and /3 = 0.65) 
is Af2oo — 5.4^^ X 10^'' Mq, again in good agreement 
with the lensing result Afaoo = ^-^-^l x 10" ^o- 

3.3.19. XLSS 0223-0436 (z = 1.22; 

XLSS 0223-0436 was discovered in the XMM Large 
Scale Structure (LSS) survey (Pierre et al. 2004; An- 
dreon et al. 2005), and Bremer et al. (2006) presented 
detailed analysis of the X-ray and optical/near-IR data 
of the cluster. They estimated an X-ray temperature of 
3.8 keV with a 1 — cr lower limit of 1.9 keV and an un- 
constrained upper limit. This X-ray temperature 3.8 keV 
is somewhat lower than what one predicts from the dy- 
namical velocity dispersion 799 ± 129 km s~^ (based on 
24 rcdshfits). Because Bremer et al. (2006) report that 
the existing XMM-Newton data do not constrain the up- 
per limit of their temperature, this cluster is excluded in 
our investigation of the lensing mass-temperature rela- 
tion. Our lensing analysis predicts — IOIII79 km s~^, 
which is even higher than the dynamical measurement. 
We estimate that the virial mass of the cluster is A/200 = 
7.4t^-5 X 10^4 Mq. 

Bremer et al. (2006) reported that the optical image of 
the cluster shows a compact distribution of > 12 galaxies 
within a 125 kpc radius centered on the X-ray peak. The 
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location of our weak-lensing mass peak agrees well with 
the X-ray and the optical centers. 

3.3.20. RDCS J0849+4452 (z = 1.26; 

RDCS J0849+4452 was discovered in the ROSAT Deep 
Cluster Survey (RDCS; Rosati et al. 1998), and fol- 
lowup near-IR imaging showed an excess of red (1.8 
< J — K < 2.1) galaxies around the peak of the X-ray 
emission. Cluster galaxies were spectroscopically con- 
firmed using the Keck telescope (Rosati et al. 1999). 
Stanford et al. (2001) measured an X-ray temperature 
Tx = 5.8^^'y keV from the Chandra data analysis; based 
on a newer calibration data, Ettori et al. (2009) reported 
Tx = 5.2 ± 1.3 keV. The cluster was one of the ACS 
GTO high-redshift cluster targets, and the color mag- 
nitude relation and the weak-lensing analysis have been 
presented by Mei et al. (2006a) and Jee et al. (2006), 
respectively. Jee et al. (2006) quoted a projected mass 
of Mproj = (2.1 ± 0.7) X 10" Mq within a 0.5 Mpc aper- 
ture radius, which is consistent with the projected X-ray 
mass ~ 2.3 x 10^"* M©, which was obtained from their re- 
analysis of the Chandra data. In the current paper, we 
perform a weak-lensing analysis using a slightly deeper 
set of the ACS images, the result of stacking both the 
GTO data (Prop. ID 9919) and the high-z supernova 
search data (Prop. ID 10496). The new mass is in good 
agreement with the Jee et al. (2006) result; in the current 



Rosati et al. (2004) presented the first X-ray analysis of 
the cluster based on both deep Chandra and XMM — 
Newton observations, which gives an X-ray temperature 



paper we quote a virial mass M; 



200 = 4.4+^;^ X 10" Mq 
(''200 = 0.98) rather than a projected mass. The cluster 
shows a strong-lensing features around the BCG. The 
two-dimensional mass map (Figure 22 ) reveals two sig- 
nificant mass clumps. The stronger one coincides with 
the centers of the X-ray emission and the optical cen- 
ter while the weaker seems to be associated with a fore- 
ground group. Both clumps were also shown in the mass 
reconstruction of Jee et al. (2006). 

3.3.21. RDCS J0910+5422 (z = l.li; 

Discovered in the RDCS, the cluster was confirmed 
with near-IR and spectroscopic observations by Stan- 
ford et al. (2002). Ettori et al. (2009) measured 
an X-ray temperature of 6.4 ± 1.4 keV from 200 ks 
Chandra data, slightly lower than the measurement 
{7.2t\l keV) by Stanford et al. (2002). Combining 
6.4 ± i.4 keV with = 147 kpc and /? = 0.843, we ob- 
tain M200 = 7.7t^:2 X 10" Mq assuming an isothermal j3 
profile. Our weak-lensing analysis yields a cluster mass 
M200 = S.Oii:^ X 10" Mq (raoo = 1-07 Mpc), slightly 
lower than, but statistically consistent with this X-ray re- 
sult. The velocity dispersion estimated from 25 redshifts 
by Mei et al. (2006) is 675 ± 190 km s-\ lower than 
our lensing prediction of 895^5^ km s^^ and the X-ray 
mass. Our mass reconstruction reveals two significant 
peaks. The stronger one is in good spatial agreement 
with the optical and X-ray centers whereas no apparent 
red-sequence is found at the location of the secondary 
peak. 

3.3.22. RDCS J 1252-2927 (z = 1.24,) 

RDCS J1252-2927 was confirmed as a cluster at z = 
1.24 based on an extensive spectroscopic campaign using 
the Very Large Telescope (VLT) (Lidman et al. 2004). 



of Tx = 6.0^0 5 keV. A revised temperature of Tx = 
6.5 ± 0.5 keV is reported in Lombardi et al. (2005) after 
the application of the new calibration of the Chandra 
instrument. The most recent analysis by Ettori et al. 
(2009) gives Tx = 7.6 ± 1.2 keV. 

Assuming the structural parameters rc = 79 kpc and 
(3 = 0.53, we translate the temperature = 7.6 ± 1.2 
keV into a virial mass M200 = (4.4 ± 1.0) x lO" Mq or a 



projected mass at r = 1 Mpc Mp 



7.1±1.1X10"M(;) 



Using ACS data, Lombardi et al. (2005) quote a pro- 
jected weak-lensing mass of M(< IMpc) = (8.0 ± 1.3) x 
10^^ Mq, which is higher than the X-ray result by ~ la. 
Our lensing analysis based on deeper ACS images gives 
a virial mass of Af2oo = 6.81|^ q x 10^^ Mq or a pro- 
jected mass of M{< IMpc) = (8.4 ± 1.2) x lO" Mq, 
which is consistent with the result of Lombardi et al. 
(2005). The current lensing analysis predicts cr„ = 
957^4g km s^^, higher than the dynamical measurement 
cr^ ^ 747ig4 km s^^ that Demarco et al. (2007) ob- 
tained from 38 redshifts. Lombardi et al. (2005) dis- 
cussed the possibility of a systematic overestimation of 
the mass in lensing because of a line-of-sight contami- 
nation. Demarco et al. (2007) reported that there is 
a possible group centered at z = 0.74 composed of 33 
members in the range 0.70 < z < 0.79. The projected 
distribution, however, is not compact, but is mostly scat- 
tered across the field. In addition, the deep 190 ks) 
Chandra data do not hint at the presence of any signif- 
icant intervening structure. 

Our mass reconstruction reveals three significant mass 
peaks. The strongest one coincides with the optical and 
the X-ray centers. The other two weaker clumps do not 
correlate well with the red-sequence of the cluster nor 
the group at 2; = 0.74. We note that these two mass 
clumps are also clearly seen in Figure 9 of Lombardi et 
al. (2005). 

4. COMPARISON WITH OTHER CLUSTER PROPERTIES 

The weak-lensing masses presented in the current work 
are used to study mass versus observable relations in the 
high-redshift universe. A comparable study but for a 
lower redshift sample was carried out by Hoekstra (2007) 
for 20 X-ray luminous clusters. Our study enables us to 
examine if the mass versus mass-observable relation holds 
across a wide range of redshifts. 

As many clusters in our sample are the results of rel- 
atively recent discoveries, a large fraction of them still 
lack information on their velocity dispersion and X-ray 
properties. In the subsequent analysis, we enlarge the 
current sample by adding five more high-redshift clus- 
ters in our previous lensing studies: XMMU J2235.3- 
2557 z = 1.4 (Jee et al. 2009), CL J1226+3332 z = 0.89 
(Jee & Tyson 2009), CL J0152-1357 z = 0.84 (Jee et al. 
2005a), MS 1054-0321 z = 0.83 (Jee et al. 2005b), and 
RX J0848-^4453 z = 1.27 (Jee et al. 2006). We refer 
readers to individual papers for detailed information for 
each cluster. 

In order to avoid any potential scatter introduced 
by variations in the mass determination method, we 
re-calculated the NEW parameters for these additional 



Scaling Relations and Overabundance of Massive High-Redshift Clusters 



9 



clusters using the same mass-concentration relation em- 
ployed in the current paper. 

4.1. Dynamical Velocity Dispersion 

Table 2 lists the dynamical velocity dispersions and the 
predicted values from the lensing analysis for the com- 
bined sample. There are a total of 23 clusters whose ve- 
locity dispersions are available either in the literature or 
through our collaborations (Meyers et al. 2011; Gilbank 
et al. in prep.). 

Figure [25] shows the comparison of the lensing predic- 
tion from the SIS fit with the dynamical measurement 
for these 23 clusters. The filled and open circles repre- 
sent the clusters analyzed in the current paper and in 
our previous studies, respectively. The dotted line is a 
fit to the data, whose slope a = 1.12 ± 0.31 and inter- 
cept 6 = 28 ± 260 km s"^ are consistent with the line of 
equality (solid); the shadow represents the 1-a range of 
the slope. Individual data points have rather large scat- 
ters around these lines mostly due to their large statisti- 
cal (e.g., small number of known spectroscopic members 
and source galaxies) and systematic (e.g., cluster mass 
profile being different from SIS) uncertainties. 

4.2. X-ray Properties 

Scaling relations between X-ray properties and lens- 
ing masses not only provide invaluable insight into the 
physical mechanism of cluster formation, but also help 
us to calibrate mass estimates based on X-ray data. For 
a virialized cluster, we expect the mass to scale with the 
gas temperature via the power law relation M cx T^/^ 
(Kaiser 1986). In this paper, we use the following spe- 
cific form to fit the data: 



E{z)Ma = A'h 



T 



5 keV 



(7) 



where E{z) is the redshift-dependent Hubble parameter 



H{z) 
Ho 



(8) 



In Equation |7j the mass Ma is usually defined as the 
total mass within the radius, at which the mean den- 
sity becomes A times the critical density. We quote 
the results for M2500 for an easy comparison with previ- 
ous studies. In our combined sample, 14 clusters have 
published X-ray temperatures (Table 2) spanning the 
1.7 keV < T < 10.4 keV range with reasonable con- 
straints on their uncertainties. For these clusters, we 
present the mass-temperature relation in Figure 26 The 
thick red line shows the best-fit power-law relation with 
a = 1.54 ±0.23 and M5 = (9.13 ±0.85) x 10" h'^ Mq. 
The cluster RCS 0439-2904 (marked with red circle) is 
a significant outlier from this relation and the fit shown 
here is obtained without including this cluster. The ex- 
clusion of this cluster in our estimation of the M-T rela- 
tion is justified because there is a strong indication that 
RCS 0439-2904 mi ght be a line-of-sight superposition of 
multiple clusters (j |3.3[ ). However, we stress that adding 
RCS 0439-2904 does not change our results (reducing the 
slope to 1.48 ±0.23). 

Allen et al. (2001) studied six relatively relaxed clus- 
ters at z < 0.46 observed with Chandra and obtained 



the mass-temperature relation a — 1.51 ±0.27 and A/5 = 
(1.32 ±0.13) X 10^4 /i-i Mq. Arnaud et al. (2005) used 
six nearby {z < 0.15) T > 3.5 keV clusters observed with 
XMM-Newton, and measured the slope a — 1.51 ± 0.11 
and the normalization M5 = (1.25±0.04) x 10" h'^ Mq, 
which are consistent with the results of Allen et al. 
(2001). A larger sample from Chandra has been studied 
by Vikhlinin et al. (2006), who analyzed 13 low-redshift 
clusters in the temperature range 0.7 — 9 keV. They also 
agreed that the observed scaling relation is in good accor- 
dance with the theoretical prediction, albeit with slightly 
higher slope 1.64 ± 0.06 than previous X-ray results. 
Their normalization M5 = (1.25 ± 0.05) x lO^ft.-^ Mq 
is in good agreement with previous X-ray results. 

Because X-ray studies derive cluster masses from tem- 
perature with hydrostatic equilibrium assumption, these 
results showing that the correlation is tight and the 
power law slope of the M — Tx relation is close to the 
theoretical prediction 3/2 may not be a direct proof for 
the cluster self-similarity. It is critical to investigate if 
the relation still holds when the masses are given by an 
independent estimator. 

The weak-lensing mass versus X-ray temperature re- 
lation is investigated by Hoekstra (2007), who derived 
a = 1.34l°|0 and M5 = (1.4±0.2) x lO^^^ h~^ Mq from 
the 17 clusters at 0.17 < z < 0.54, which are in good 
agreements with the values from X-ray studies. Mahdavi 
et al. (2008) updated the redshift distribution used by 
Hoekstra (2007), and found that the new n{z) from the 
photometric redshift catalog of Canada-France Hawaii 
Telescope Legacy Survey (CFHTLS) derived by Ilbert et 
al. (2006) gives on average ^10% smaller masses. We 
show this revised result in Figure [26] We further note 
that Hoekstra (2007) used the ASCTv temperatures that 
do not correct for cool cores. The absence of this cool- 
core correction can bias the temperatures low. 

We verified that the above M — Tx relation is not sen- 
sitive to a mass-concentration relation. When we recom- 
pute our A/2500 using the mass-concentration relation of 
Bullock et al. (2001), the power-law slope a virtually 
remains the same and the normalization (M5) decreases 
only ~3%. The X-ray and Hoekstra (2007) results were 
obtained without any assumption on mass-concentration 
relation. Nevertheless, we tested if the result of Hoek- 
stra (2007) changes when the masses are derived from 
NFW profile fitting with the Duffy et al. (2008) mass- 
concentration relation as is done in this paper. Again, 
we found that neither the normalization nor the slope of 
the M — Tx relation of Hoekstra (2007) changes because 
of this mass determination met hod. 

It is easy to see in Figure [26] that there is a 20 — 30% 
discrepancy in normalization between our result and pre- 
vious results. The discrepancy may imply that the nor- 
malization decreases with redshift or our weak-lensing 
masses are biased low. One cause for possible weak- 
lensing bias is the redshift estimation bias, similar to 
the case in Hoekstra (2007). We estimate that the sam- 
ple variance of UDF may be responsible for 5-12% shift 
in mass depending on the cluster redshift. However, al- 
though we cannot exclude this possibility, this would ag- 
gravate the already problematic existence of the most 
massive clusters in the high-z universe (see y5]for details). 
Therefore, it seems more plausible that tnere is some 
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evolution in the normalization of the mass-temperature 
relation to the highest redshifts. 

Finally, we compare the masses derived from X-rays 
and lensing analyses. The cluster masses at r2oo obtained 
from the X-ray results are listed in Table 2. We use an 
isothermal j3 assumption to derive these values in order 
to take advantage of the published parameters in the 
literature. For the clusters that do not have published 
measurements of their X-ray surface brightness profile 
(because the S/N of the existing data do not constrain 
the shape of the X-ray profile well) we assume /3 = 0.7. 
Figure [27] shows the result for 14 clusters in our com- 
bined sample. The agreement between the lensing and 
X-ray estimates is good except for RCS0439-2904, which, 
as mentioned already, might be a line-of-sight superpo- 
sition of multiple components (Cain et al. 2008). It is 
important to remember that an SIS profile is assumed 
for the derivation of M200 from X-ray data. Therefore, 
it is possible that this excellent agreement might be a 
coincidence to some extent and a systematic difference 
might appear if masses are evaluated at different radii. 

5. ABUNDANCE OF THE MOST MASSIVE CLUSTERS AT 
HIGH REDSHIFT 

The abundance of the most massive clusters at high 
redshift is extremely sensitive to cosmological parame- 
ters. In the current era when we can estimate the cos- 
mological parameters to a high precision [mainly from 
the cosmic microwave background (e.g., Komatsu et al. 
2011) and Type-la supernova studies (e.g., AmanuUah 
et al. 2010)], discovery of even a single massive cluster 
can challenge the current ACDM model (Mortonson et 
al. 2011). 

In Jee et al. (2009), we presented a weak-lensing anal- 
ysis of XMM J2235-2557 and found that the mass of 
the cluster is surprisingly high. With the cosmological 
parameters fixed at the current best-fit values, the the- 
oretical probability of finding such a massive cluster in 
the survey volume is only '^1%. The exact value for 
the probability depends on the adopted mass function 
and cosmological parameters. Nevertheless, independent 
estimates (e.g., Jimenez & Verde 2009; Sartoris et al. 
2010; Holz & Perlmutter 2010) agree that the discov- 
ery provides significant tension with the current ACDM 
cosmology. In this section, we extend the study of Jee 
et al. (2009) to the most massive clusters in our sam- 
ple with a new mass function, revised survey areas, and 
marginalization over the uncertainties of cluster masses 
and cosmological parameters. 

5.1. Mass Function 
The mass function is often given in the following form: 
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where a is the rms variation of the density field when 
smoothed on scale M. Pm.o is the present matter den- 
sity. The function / determines the shape of the mass 
function. 

Jenkins et al. (2001) proposed a redshift-independent 
form / = f{o'), which makes the above mass function 
agree with their numerical simulation results with high 
precision. Tinker et al. (2008) improve the agreement 



by considering the redshift-dependence of f — f{a, z) 
which is expressed as 



/(a,z) = A 





— a 




[(^) 




hi 









(10) 



where A = 0.186(1 -h z)-^-^^, a = 1.47(1 + z)'°-°'^, 
b = 2.57(1 -I- z)"°°", and c = 1.19 when the mass is 
defined as the total mass at the radius, inside of which 
the mean density is 200 times the mean density of the 
universe at the cluster redshift. We refer to it as M2oom 
when necessary to distinguish it from M200, which is 
more commonly used to refer to the mass at the radius, 
where the mean density becomes 200 times the critical 
density. The deviation of the universality of the mass 
function by Jenkins et al. (2001) from the numerical 
simulation results is 20-50% (Tinker et al. 2008); Bhat- 
tacharya et al. (2011) argue that the difference between 
different fitting formulae is probably larger than the scat- 
ter seen in the largest simulations. 

5.2. Cluster Survey Area 

The cluster survey area that one chooses to adopt to 
evaluate the cosmological significance for discovery of un- 
usually massive clusters is somewhat controversial. For 
example, Mortonson et al. (2011) argue that ^300 deg^ 
should be used for the search area of XMM2235-2557 
at 2: = 1.4 instead of the pilot survey area 11 deg^ 
(Jee et al. 2009). The argument for this - 300 deg^ 
area is based on the study of Hoyle et al. (2011), who 
assumed -168 deg^ for XCS, ~64 deg^ for XMM-LSS, 
~11 deg2 for XDCP, -17 deg^ for WARPS, and ~33 
deg2 for RDCS. However, the different XMM-Newton- 
based surveys overlap significantly. In addition, because 
the flux limit that enables a secure cluster detection at 
z > 1 is ~ 10~^'*erg cm~^s~^, the effective search area 
must be reduced significantly. For example, the XDCS 
survey has continued since the end of the 11 deg^ pi- 
lot survey and now reaches ~ 50 deg^ (R. Fassbender in 
prep.) for a flux limit ~ lO'^Vg cm'^s'^ The RDCS 
survey covers a geometrical area of 50 deg^ only at rel- 
atively high fluxes. However the area suitable for z > 1 
cluster detection is ^5 deg^, corresponding to fluxes be- 
low ~ 3 x 10-"erg cm^^s"^ (Rosati et al. 1998). Con- 
sidering the argument listed above, the combined survey 
area for the X-ray selected high-z clusters should be ^100 
deg^. 

5.3. Abundance Estimation and Eddington Bias 

In the full sky survey, the number of clusters with mass 
and redshift greater than Mmin and Zmm, respectively is 
given by 
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where dV/dz is the volume element and dn/dM is the 
mass function. For massive high-redshift clusters, the 
result is sensitive to the mass function dn/dM near Zmm 
and Mjnin- 

We compute the probability of finding each cluster by 
marginalizing over the measurement errors, assuming a 
log-normal distribution for P{M). This is different from 
the approach of Jee et al. (2009), where we simply used 



Scaling Relations and Overabundance of Massive High-Redshift Clusters 



11 



a threshold mass Mthr to evaluate the integral (EqnfTTl). 
Nevertheless, we emphasize that the new method yields 
values very close to those obtained from the previous 
method when the mass uncertainty is relatively small 
(< 20 %). For example, marginalizing over the measure- 
ment error gives an abundance of ~ 2 x 10~^ for XMM 
J2235-2557 in the original 11 deg^ XDCP survey, which 
is in good agreement with the previous value ~ 5 x 10~^. 
As we use a fixed value of Zmax i the resulting abundance 
is slightly conservative because in principle the maximum 
redshift should decrease at the low end of P{M); remem- 
ber that at the high end of P{M) , the steep mass function 
cancels the effect of increased Zmax- For the conversion 
from expected number to probability we assume Poisso- 
nian distribution. That is, if the expected number is N 
within a particular survey, the discovery probability is 
simply P = 1 — e~^. 

Mortonson et al. (2011) noted that an Eddington bias 
(Eddington 1913) is an important factor in the discus- 
sion of the cosmological significance using most mas- 
sive clusters. Eddington bias is the combined effect 
of a steep mass function and measurement uncertainty. 
For example, consider the weak-lensing mass estimate of 
XMM2235 M200 = (7.3 ± 1.3) x IO^Mq. Because the 
slope of the mass function near the cluster mass is steep, 
there should be more < 7.3 x IQ^^Mq objects scattering 
up than > 7.3 x IO^^Mq objects scattering down. As a 
matter of course, this does not mean that XMM2235's 
mass should be quoted with a lower mass (for individual 
clusters the scatter is symmetric). However, when we 
discuss the existence of the cluster in probabilistic terms, 
this additional chance of bias should be considered from 
a Bayesian point of view. 

We adopt the convenient Mortonson et al. (2011) pre- 
scription, which is to replace the observed mass with the 
following reduced mass M' in the estimation of the abun- 
dance: M' = exp(l/27(7inM)-M', where 7 is the local 
power-law slope (rfn/dlnM ~ M^), and cinAf is the 1-cr 
uncertainty of In M (log-normal distribution is assumed 
for mass errors). 

5.4. Probability of Discovery for Most Massive Clusters 

In Figure [28] and Table 3 we show the discovery proba- 
bility of all clusters (i.e., similarly massive clusters) stud- 
ied in this paper. We compute the probabilities after 
marginalizing over the Wilkinson Microwave Anisotropy 
Probe 7-year (WMAP 7) Monte Carlo Markov Chain 
(MCMC) datc[^ The probabilities estimated within the 
parent survey may lead to a somewhat aggressive inter- 
pretation for some clusters. For example, most of the 
RDCS clusters have significantly low probabilities mainly 
due to the small (~5 deg^) survey area. However, most 
non-X-ray selected clusters (ISCS and RCS) have discov- 
ery probability equal or close to unity even within their 
parent surveys. 

As mentioned in j ]5.2[ we believe that the survey area 
100 deg^ is still a conservative choice for most of the 
clusters shown here (the exceptions are MS1054-0321 
and CL J0152-1357, which were discovered in the EMSS 
~ 735 deg^ survey). With this survey area, there re- 
main four clusters whose discoveries are still unusual 



' available at 



given our current understanding of cosmological param- 
ete rs ( the names of these clusters are shown in Fig- 
ure[28|. The cluster with the lowest probability (~1%) is 
the z ^ 0.9 cluster CLJ1226-I-3332, whose weak-lensing 
mass and X-ray temperature are M2oom ^ 1-6 x IQ^^Mq 
and Tx ~ 10 keV, respectively. Because the clus- 
ter was discovered in a 72 deg^ survey (WARPS), us- 
ing 100 deg^ does not increase the probability signifi- 
cantly. In addition, XMM-Newton archival search pro- 
grams look for regions where no known X-ray luminous 
clusters exist, and thus the search volume in the XMM- 
Newton archival survey excludes the z ^ 1 volume by 
observer's selection. This makes our probability estima- 
tion for CLJ1226-I-3332 even more conservative. The sec- 
ond lowest probability is found with XMM2235-2557 at 
z = 1.4. The probability ^3% is much higher (after 
the area and Eddington bias correction) than our previ- 
ous value ~0.5% within the 11 deg^ pilot XDCP survey, 
but still gives ~ 2a tension with the predicted number 
of clusters from ACDM. RDCS1252-2927 at z = 1.24 is 
ranked as third in rarity, giving P ^-^6%, which is followed 
by XLSS0223-0436 at z = 1.22 (P -18%). 

Combining discovery probabilities of all our clusters is 
somewhat tricky. Because of Poissonian fluctuation, the 
discovery probability is always lower than the estimated 
abundance (i.e., P = 1 — e^^). For example, if the ex- 
pected number of a given cluster is A'^ = 1 within 100 
deg^ survey, the probability is ~ 0.63. Obviously, sim- 
ply multiplying probabilities of all clusters potentially 
leads to an unreasonably low final value if there are many 
slightly-less-than-one-probability clusters. Nevertheless, 
if we include only rare clusters (P << 1), multiplying in- 
dividual probabilities provides a useful measure on how 
much the existence of these multiple clusters give rise to 
tension with the current cosmological parameters. 

One set of cosmological parameters that increase the 
predicted abundance of one massive cluster can boost an- 
other similarly massive cluster's abundance. Therefore, 
multiplication of probabilities must be done separately 
for each chain, and the final probability should be ob- 
tained after properly taking into account each chain's 
weight. We find that the final combined probability ob- 
tained in this way is —0.03%. 

5.5. Exclusion Curve Test 

Mortonson et al. (2011) presented fitting formulae, 
which provide an exclusion mass as a function of red- 
shift for given sample and parameter variance confidence 
limits (CL). Even a single cluster equal to or above the 
exclusion mass rules out both ACDM and quintessence. 

In Figure |29[ we compare our clusters to the exclusion 
curves of Mortonson et al. (2011). Although we believe 
that the proper choice for the sky area is 100 deg^ or less, 
we also plot the result for the 300 deg^ assumption, which 
was adopted in Mortonson et al. (2011). For simplicity, 
we only consider the case where the sample variance CL 
is equal to the parameter variance CL, which is referred 
to as joint CL in Mortonson et al. (2011). 

The central mass of the cluster CL J1226-I-3332 after 
the Eddington bias correction lies above the joint CL 95% 
exclusion curve when we adopt 100 deg^ for the survey 
area. All the masses of the four clusters shown here are 

http://lambda.gsfc.nasa.gov/product/map/current/parg|ig ^^igjjg|.^]tyi_Q^^ the joint CL 80% exclusion 
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curves. When we consider 300 deg^ for the search area in- 
stead, we note that the two clusters CLJ1226+3332 and 
XMM J2235-2557 are above the 80% exclusion curves. 
We clarify here that Mortonson et al. (2011) used the 
mass of XMM J2235-2557 reported by Stott et al. (2011), 
whose value is similar to our weak-lensing mass, but 
has significantly larger uncertainty 40%); Stott et al. 
(2011) derived the cluster mass from their X-ray scaling 
relation. This large uncertainty leads to ~ 40%. Edding- 
ton bias, and Mortonson et al. (2011) concluded that 
the resulting mass of the cluster provides only a negligi- 
ble tension with ACDM when they assume 300 dcg^ for 
the sky area. On the other hand, our mass imcertainty is 
~ 18%, and the corresponding Eddington bias is ^ 12%, 
which makes the corrected mass still close to the ~85% 
exclusion curve for Mortonson et al. (2011)'s 300 deg^ 
assumed area. 

The above exclusion curve test provides weaker con- 
straints when one is using more than a single cluster to 
examine their cosmological significance because the tool 
adopts the lowest redshift for the sample in the estima- 
tion of the cluster abundance (Mortonson et al. 2011). 
Therefore, we do not attempt to estimate the tension by 
combining the four most massive clusters. 

5.6. Interpretation of Overabundance of Massive 
Clusters at High Redshift 

Assuming the current well-accepted ACDM cosmol- 
ogy, we have demonstrated above that the discovery of 
the most massive clusters in our sample are rare events. 

Among the most popular interpretation would be that 
the primordial density fluctuation is non-Gaussian. For 
example, Jimenez & Verde (2009) claim that the local 
/jVL should be in the range 150-260 to accommodate 
XMMU2235-2557 within the pilot 11 deg^ survey We 
suspect that their estimate of /^vl would decrease some- 
what when both the Eddington bias and revised area 
50 — 100 deg^ are taken into account. A tighter con- 
straint is possible if one utilizes the weak-lensing masses 
of the four most massive clusters studied here. A similar 
study is performed by Hoyle et al. (2011) from the inves- 
tigation of the 14 high-redshift cluster masses reported 
in the literature, giving a significantly large f^^ > 375 
(after marginalization over WMAP5 parameters). 

One crucial test that however should accompany these 
non-Gaussianity studies is the investigation of the mass 
function at the high mass end. The existing limited- 
volume simulations do not constrain the number density 
of extremely massive clusters accurately, and the com- 
monly used fitting functions are simply extrapolated re- 
sults in this regime. Another potentially important con- 
tribution from numerical simulation is the predicted mass 
function specific for a given survey, which takes into ac- 
count the various aspects of selection limits and projec- 
tion effects. 

The projection effect is always a concern in the cosmo- 
logical interpretation of extremely massive clusters. Cer- 
tainly, a superposition of two moderately massive clus- 
ters or a long filament viewed along the line of sight 
can be identified with an extreme object. The former 
case is easily detected by a redshift histogram if the field 
is spectroscopically surveyed. Initially, RDCS1252-2927 
at z = 1.24 was such a candidate because the redshift 
catalog showed three possible foreground groups at at 



z = 0.47, 0.68, and 0.74 (Lombardi et al. 2005). How- 
ever, the spatial distributions of these groups are not 
compact, and they appear to be just loose galaxy groups 
commonly found in any random field. The latter case (a 
filament viewed along the line of sight) is very difficult 
to prove (or disprove) directly. Nevertheless, there are 
a few circumstantial lines of evidence, which help us to 
infer if a system possesses extreme triaxiality. Compar- 
ing X-ray temperature with weak-lensing mass is use- 
ful because a severe triaxiality will give rise to much 
higher weak-lensing mass than what the temperature 
predicts. In addition, a line-of-sight elongation inflates 
velocity dispersion, and the result will be significantly 
larger than what the Icnsing predicts. XLSS0223-0436 
is potentially such an object because our lensing mass 
M200 = 7.4ti | X 10^^ Mq is much larger than what the 
X-ray temperature ^ 3.8 kcV (Bremer et al. 2006) sug- 
gests. However, the X-ray photon statistics are poor, and 
Bremer ct al. (2006) could not determine an upper limit 
of the X-ray temperature. Deeper Chandra observations 
would be necessary to robustly examine this large (a fac- 
tor of 2-3) discrepancy. We note that the existing data 
for CLJ1226+3332 and XMM J2235-2557, the two clus- 
ters with the largest tension with ACDM in our sample, 
do not suggest the possibility that weak-lensing masses 
are significantly boosted by any projection effect. 

6. SUMMARY AND CONCLUSIONS 

We have presented a weak-lensing analysis of 22 high- 
redshift (0.91 < z < 1.46) clusters using' HST/ ACS im- 
ages. Most clusters in the sample show a significant lens- 
ing signal and are well detected in our two-dimensional 
mass reconstruction. We complement our sample with 
five high-redshift clusters that were studied in the past 
in order to compare lensing results with X-ray and dy- 
namical measurements over a large range of masses, and 
to discuss the most massive clusters in this combined 
sample in the context of cosmology. Our results are sum- 
marized as follows. 

1. The dynamical velocity dispersions are consistent 
with the values estimated by lensing analysis. 

2. The power law slope of the lensing mass- 
temperature relation is in good agreement with the 
theoretical prediction, which is also supported by 
previous investigations carried out with X-ray data 
alone or with both X-ray and lensing data. 

3. We see evidence of an evolving normalization con- 
stant in the M — Tx relation. 

4. The masses derived from X-ray and lensing mea- 
surements are in good agreement. 

5. The existence of the most massive clusters in our 
sample gives rise to tension with the current ACDM 
standard structure formation paradigm. 

The first conclusion is consistent with the findings of 
Hockstra (2007) who studied low redshift clusters. Our 
result demonstrates that the relation can be extended to 
this high-redshift regime, albeit with slightly increased 
scatter. This larger scatter could be intrinsic for high 
redshift clusters, which are likely to possess more merging 
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substructures than low redshift systems. However, one 
must remember that most clusters in our sample do not 
have sufficient number of spectroscopic members yet to 
give reliable estimates. In addition, the depth of the ACS 
images is also less than optimal in many cases, which in- 
creases uncertainties for the lensing values. Our assump- 
tion of SIS in the derivation of the velocity dispersion 
may introduce some bias in the relation simply because 
the real clusters are not singular isothermal. Neverthe- 
less, the collective behavior of the relation being consis- 
tent with the equality is still noteworthy, and hints at 
the possibility that the dynamical measurements from a 
large number of cluster members may still serve as usehil 
proxies for masses in this redshift regime. 

The mass-temperature relation presented here is tight, 
particularly for hot clusters. The current result indi- 
cates that the power-law slope of the mass-temperature 
relation does not evolve since z ~ 1. However, our nor- 
malization value is 20 — 30% lower than previous results, 
suggesting that the X-ray temperatures of the clusters at 
z ~ 1 may be systematically higher than at low redshifts. 
The shift can imply that non-thermal sources of heating 
(e.g., merging, AGN activity, etc.) might be more preva- 
lent at higher redshift, or simply the difference in gas 
fraction (because of different cluster formation epoch) 
might be significant. Alternatively, one can also consider 
bias in our lensing analysis. However, it is very unlikely 
that our calibration error, even at a maximum, can in- 
duce 20 — 30% systematic errors. Our shear calibration 
error verified by numerical simulations is at most at the 
^ 1% level. The redshift estimation bias (i.e., the magni- 
tude versus redshift relation in the UDF does not repre- 
sent the mean relation of the universe) is at the 5 — 12% 
level for each cluster, giving the highest error for the most 
distant object. Most importantly, if the cluster masses 
in our sample should be somehow increased by 20 — 30%, 
this would not only degrade the dynamic versus lensing 
velocity dispersion relation substantially, but also make 
the existence of the most massive clusters much more 
problematic with the current cosmological parameters. 

The fourth conclusion is somewhat foreseen by the 
tight mass-temperature relation. However, this direct 
comparison of mass between X-ray and lensing methods 
provides another way to verify our normalization in the 
M — T relation. Again, the clusters with the lowest tem- 
peratures are the greatest outliers from the equality. On 
the other hand, the other relatively massive clusters give 
consistent results between X-ray and lensing measure- 
ments. When we repeat the analysis by re-calculating 
the cluster mass at smaller radii (e.g., M500 instead of 
-M200), the relation virtually remains the same. This in- 
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dicates that the difference in the assumption of the clus- 
ter mass profile (i.e, isothermal [3 versus NFW) is not a 
significant source of bias between the two measurements. 

Finally, the last conclusion is the most striking find- 
ing of the current investigation. The above studies on 
the relations between lensing mass and other properties 
do not support the possibility that our lensing mass de- 
termination is significantly biased high. A significantly 
high degree of non-Gaussianity or high value of (Js may 
explain the existence of these massive clusters at high 
redshift. However, we caution that our current theoreti- 
cal knowledge of the mass function at high end is limited 
by small number statistics. 

Foley et al. (2011) reported the discovery of SPT- 
CL J2106-5844, which is claimed to be the most massive 
cluster known at 2; > 1 based on the X-ray temperature 
and SZ measurement. The authors estimate that there is 
only a 7% chance of finding such a massive cluster within 
the 2,500 dcg^ survey area. The mass of the cluster still 
needs to be verified by other methods (e.g., weak lens- 
ing analysis). However, if the cluster is indeed massive 
as reported, the existence of SPT-CL J2106-5844 also 
favors our claim that the most massive clusters at high 
redshift provide non-negligible tension with our current 
understanding of structure formation. 
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Fig. 1. — Reduced tangential shears for 22 high-redshift clusters. Filled and open circles represent the tangential shear and 45 ° rotation 
test results, respectively. 
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Fig. 2. — Reduced Tangential shears - continued from Figure [T] 
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Fig. 3. — XMMXCS J2215-1738 mass map. The left panel displays the mass reconstruction, whose color represents the mass density. 
The mass-sheet degeneracy is not broken here, and thus the scale is arbitrary. On the right panel, we overlay the mass contours on the 
pseudo-color composite created by combining the 1775 and zgso images. The green "X" symbol marks the location of the X-ray peak in the 
Chandra image. Spectroscopically confirmed members are encircled in red. 




Fig. 4.— 



Same as Figure|3]but for XMMU J2205-0159. 




Fig. 6.— 



Same as Figurc|3]but for WARPS J1415+3612. 




Fig. 7.— Same as Figure|3]but for ISCS J1432+3332. 
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Fig. 9.— 



Same as Figurc|3]but for ISCS J1434+3427. 
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ame as Figure [s] but for ISCS J1434+3519. 




Fig. 13.— 



Same as Figure [s] but for RCS 0220-0333. 
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ame as Figure [s] but for RCS 0337-2844. 




Fig. 17.— 



Same as Figure [s] but for RCS 2156-0448. 





ame as Figure [s] but for RCS 2345-3632. 
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Same as Figure [s] but for XLSS J0223-0436. 



Scaling Relations and Overabundance of Massive High-Redshift Clusters 



O.™ aOJi 0.07 0.10 0.13 D-.I7 0.39 




ame as Figure |3] but for RDCS J0849+4452. 
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ame as Figure [s] but for RDCS J0910+5422. 
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Fig. 24.— Same as Figure |3] but for RDCS J1252-2927. 
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Fig. 25. — Dynamical velocity dispersion versus lensing prediction. We use filled circles to represent the results for the clusters studied 
here. Open circles show the results for the high-z clusters in our previous publications. The solid line represents the line of equality The 
lineaj: fit to the result is shown with the dotted line while the 1-a range of the slope is depicted by the grey region. 
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Fig. 26. — X-ray temperature versus lensing mass. As in Figure pS] filled and open circles represent the current clusters studied here 
and the five clusters in our previous publications. The slope of the power law ct = 1.54 it 0.23 (M oc T") is consistent with the theoretical 
prediction 3/2 and also previous results obtained from low redshift samples. Note that the Hoekstra (2007) relation shown here is the 
revised result after we applied the 10% reduction in mass as explained in Mahdavi et al. (2008). The greatest outlier from the M — Tx 
relation is RCS 0439-2904 (red circle), which perhaps is a line-of-sight superposition of multiple components (Cain et al. 2008) boosting 
the weak-lensing mass measurement without significantly affecting temperature measurement. 
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Fig. 27. — Virial mass comparison betw een X-ray and weak-lensing results. Not surprisingly, the trends are similar to that of the mass- 
temperature relation shown in Figure [26] The cluster RCS0439-2904 (red circle) again is the greatest departure from the line of equality. 
Filled and open circles represent the current clusters studied here and the five clusters in our previous publications, respectively. 
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Fig. 28. — Discovery probability of clusters. After applying Eddington bias to the lensing masses, we marginalize over the uncertainties 
of the WMAP7 cosmological parameters in order to estimate the discovery probability. MS 1054-0321 and CL J0152-1357 were discovered 
in the Einstein Medium Sensitivity Survey (EMSS), which covered a total of 778 deg^ (Gioia et al. 1990). We do not recompute the 
probability for these two clusters with 100 deg^. 
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Fig. 29. — Exclusion curves. Any cluster's mass above the given curve excludes the ACDM with the joint CL specified (Mortonson et al. 
2011). For example, the central value (Eddington bias-corrected) of CL J1226-I-3332 rules out >95% of currently allowed ACDM models 
at the > 95% confidence level for the 100 deg^ area. 
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Details of HST observations 
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Cluster Properties 
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-0.13 
+0.12 
-0.11 
+0.77 
-0.07 
+0.08 
-0.07 
+0.06 
-0.06 
+0.08 
-0.07 
+0.10 
-0.09 
+0.11 
-0.09 
+0.09 
-0.06 
+0.09 
-0.07 



4 3+3° 



3.0: 

5.3 
4.7- 
4.9 
5.4 
2.5 
5.3 
2.8 



+1.6 
1.0 

+1.7 
1.2 

+2.0 
-1.4 
+ 1.6 
1.2 
+2.4 
1.6 
+2.2 
1.1 
+2.6 
1.7 
+2.9 
1.4 



3.1^ 



1+2.6 
-1.4 

4 8+1-* 

1 0+1.3 
1-8-0.7 

4.9tli 
A 3+1-7 
^-•^-1.2 
1 0+2.5 
1-8-1.0 

1 9+1'' 

2 4+11 

5 8+2-3 

7 4+2-5 
'-^-1.8 
4 4+1-1 
^■^-0.9 

5- 011:^ 

6- 811:^ 
7 0+1.7 

+2.4 
2.0 
+2.1 
1.8 
4 4+0-7 

+1.0 
0.8 



13.7 
10.8 



3.1"' 



Note. — " Hilton et al. (2010). '' Santos et al. 
al. (2011). 9 Andreon et al. (2008). Cain ct al. 
" Demarco ct al. (2007). " Rosati ct al. (2009). 
Santos et al. (2009). ^ Maughan et al. (2006). ' 

Ettori et al. (2009). ® Rosati et al. (2009). Maughan et al. (2007) 



(2009). ■= Huang et al. 

(2008). • Gilbank ct al. 
° Holdcn ct al. (2009). 
Andreon et al. (2011^ 



(2009). 



Eisenhardt et al. (2008). ' Meyers et al. (2011). f Brod-win et 

in prop. ^ Gilbank ct al. (2008). * .Icc ct al. (2006). ' Moi ct al. (2006). 
P Tran ct al. (2007). ' Demarco ot al. (2005). ^ Hilton ct al. (2010). ^ 
^ Cain et al. (2008). ® Hicks et al. (2008). ^ Bremer et al. (2006). 
Jee et al. (2005b). 
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TABLE 3 

Discovery Probability of Galaxy Clusters 



Cluster name 


Within Parent Survey 


100 deg2 


AlVlMACb JZzlo-iYoo 


0.96 


1 


aMMU JzzUo-Uloy 


1 


1 


aMMU Jlzzy+Uiol 


0.61 


0.81 


WAHFb J 1415+0012 


0.65 


0.96 


TQ/~iQ 11 /I QO 1 QQQO 

IbUo J14oi+oooz 


0.14 


U. / / 


T O C T"1 /ion 1 O/ID"? 


0.15 


0.79 


TC/~^C T "1 A A 1 O /I T7 

ibCb J14d4+J427 


1 


1 


T O C T1 /I OO 1 O /I O/^ 

ibCb J14di+o4do 


0.11 


0.70 


IbOb J 1434+3519 


1 


1 


IbL/b J14ao+o414 


0.92 


1 


T>/~<0 r\oor\ noQO 


0.74 


0.91 


T>/~iO r\ooi r»ooi 

KCb Uzzl-Uozl 


1 


1 


HOb Uc>c)7-io44 


0.84 


1 


T~>/^o n /I on nnn /I 

HOb U4c)y-iyu4 


0.95 


1 


T> /^o nice n/i/io 

KCb Z100-U448 


1 


1 


rtOo loilH-uyUo 


1 
i 




RCS 2345-3632 


1 


1 


RCS 2319+0038 


0.83 


0.93 


XLSS J0223-0436 


0.01 


0.18 


RDCS J0849+4452 


0.03 


0.60 


RDCS J0910+5422 


0.06 


0.65 


RDCS J1252-2927 


0.002 


0.052 


XMMU J2235-2557 


0.013 


0.024 


CL J1226+3332 


0.006 


0.01 


MS 1054-0321 


0.35 


_a 


CL J0152-1357 


1 


_a 


RDCS J0848+4453 


0.08 


0.88 



Note. — " MS 1054-0321 and CL J0152-1357 were discovered in the Einstein Medium Sensitivity Survey (EMSS), which covered a total of 778 
deg^ (Gioia et al. 1990). 
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APPENDIX 
A. PSF MODELING 

The importance of accurate removal of PSF effects cannot be stressed too much in weak-lensing analyses. Both 
dilution and bias of shear signals must be accounted for with rigor in order to advance from mere detection into 
measurement. One difficulty in ACS PSF modeling is the small field of view (3' x 3'), which does not provide a sufficient 
number of high S/N stars. Because ACS PSFs vary with time and position, it is not possible to interpolate/extrapolate 
the PSF information reliably based a few tens of stars on a target image. Fortunately, the PSF pattern in ACS seems 
to be repeatable. That is, one can find two or more observations remote in time but closely related to each other in 
the behavior of the PSF pattern (Jee et al. 2007b). Therefore, it is possible to use a limited number of stars on one 
image, find another image possessing a similar PSF pattern but instead with many more stars covering the entire field, 
and utilize this second image to infer the PSF in the first image. In practice, a slight complexity arises because in 
general a target field is visited with many telescope pointings that are different not only in shift and rotation, but also 
in PSF pattern. Therefore, one needs to determine PSF patterns for individual pointings first, and then apply shift 
and rotation to model the PSFs on the final mosaic image. 

The above method has been successfully applied to our previous HST weak-lensing analysis of clusters. We refer 
readers to those papers (e.g., Jee et al. 2009) for details. Here, we only summarize the PSF modeling procedure 
for the current data sets. We use the Jee et al. (2007b) PSF library, which provides a polynomial description of 
two-dimensional PSF variations based on principal components analysis (PC A). This PCA method allows us to obtain 
compact basis functions to describe the spatial variation of PSFs (Jee et al. 2007; Jee & Tyson 2011). 

We identified stars on each cluster observation using objects' size and magnitude. The best PSF template for each 
cluster observation frame was determined by comparing both size and ellipticity of these stars. We stored the name 
of the template, and the alignment information (i.e., shift in pixels and rotation angle in degrees). After completing 
this task, we then went through objects in the source catalog that is created from the final stack, computed their 
location within each input frame, and interpolated the coefficients of the principal components at the location. These 
coefficients were transformed to two-dimensional PSF images, which were stacked to generate the final PSF model 
after proper rotation was applied; in fact, the charge transfer inefficiency (CTI) correction (see j |B| was made prior to 
stacking by stretching the individual PSFs in the direction of CTI trail. 

In Figure [30]|35] we show the PSF correction results. We note that the observed PSF ellipticity on the stack is 
on average smaller than the typical value because of the rounding effect from field rotation. The rms residual per 
cluster is less than 1%. Despite some occasional small mismatches, the current level of the PSF correction exceeds the 
requirement for cluster lensing analysis. 

B. CTI CORRECTION 

During a CCD readout, some fraction of the charge is trapped because of defects in the silicon. These trapped 
charges are soon released after a characteristic time r. This cycle of trap-release events continues until the charges in 
the last row are transferred. This artifact is sometimes clearly visible as "charge trails" . 

The HST/ACS CTI worsens with time because the number of defects in the CCD is proportional to the amount 
of the time exposed to the radiation-rich space environment. Most of the cluster images in our program are obtained 
during the year 2006 (four years after the installation of ACS), and thus we expect the CTI to be much worse than in 
our previous studies. 

Obviously, the effect is undesirable in most applications of the data. The effect is particularly undesirable in weak- 
lensing measurements, where one is looking for a subtle distortion of small galaxies. In this analysis, the CTI-induced 
elongation must be carefully measured and accounted for. 

In the characterization of the CTI, we categorize the methods in the literature into four schemes. 

1. Difference photometry. 

2. Pixel Response measurement. 

3. Ellipticity bias of astronomical objects. 

4. Ellipticity bias of sub-seeing features (e.g., cosmic ray, hot or warm pixels, etc.). 

5. Pixel level correction. 

Method 1 utilizes the fact that the CTI-induced elongation leaves some fraction of the charge outside the aperture. 
Because the fractional loss depends on the CTI, in this scheme an object should be observed multiple times in 
different conditions (e.g., background, time, row, column, etc.). Therefore, this method requires large amounts of well- 
planned, dedicated observations (e.g., Riess & Mack 2004). This scheme is also useful when one's interest is limited to 
photometric correction of the CTI effect. Method 2 involves uniformly illuminating a CCD and measuring deviations 
from that uniformity in the few pixels farthest from the readout register. A series of programs for First Pixel Response 
(FPR) and Extended Pixel Edge Response (EPER) monitoring have been carried out (Prop. ID 9649, 10044, 10369 
and 10732). A similar, but more elaborate measurement of this kind is possible if the CCD is exposed to a number of 
point sources (i.e., single pixel events created by a controlled X-ray source). However, this radiation test is only possible 
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on the ground. Method 3 and 4 rely on the fact that for ACS, the CTI-induced elongation is mostly significant along 
the parallel charge-transfer direction (i.e., y-axis in CCD coordinates). Therefore, taking the average over a number 
of measurements should reveal the net CTI effect. When employing Method 3, one has to ensure that the shapes of 
the employed objects are free from other systematic effects (e.g., Schrabback et al. 2010; Hoekstra et al. 2011). For 
example, poor PSF modeling can induce a residual shape bias that is similar to the CTI-induced distortion in both 
direction and flux-dependence. Method 4 differs from Method 3 in that the ellipticity bias is measured from cosmic 
rays or warm pixels, whose shapes are not subject to geometric distortion or PSF. Because these sub-seeing features 
are numerous even on a single 500 s exposure image, it is possible to obtain high S/N information on CTI from the 
data themselves. We applied this method to the weak-lensing analysis of the then highest redshift cluster XMM2235 
at ^ = 1.4 (Jee et al. 2009). As Jee et al. (2009) were somewhat implicit on how to relate the measured ellipticity 
bias of cosmic rays to object ellipticity, we provide details in the current paper. Method 5 is recently suggested by 
Massey et al. (2010) and Anderson & Bedin(2010). This approach aims to restore CTI-trails back to the original pixels 
based on profile measurements of warm-pixels. Potentially, this method provides convenient CTI-corrected images, 
which astronomers can work on directly without worrying about the details of the model. However, The fidelity of the 
method has yet to be tested. 

B.l Quantification of CTI Effect 

Any significant sub-seeing features (hereafter SSFs) present in the raw CCD images are mostly cosmic rays, warm/hot 
pixels, or simply photon noise, which are commonly removed through elaborate image-processing. However, knowing 
that these SSFs are not caused by the photons following the regular optical path of the instrument, we can use these 
otherwise nuisance features to our advantage. Because position angles of the SSFs in principle should possess no 
preferred direction on the surface of the CCD, any deviation from the anticipated isotropy is a clean indicator of the 
artifact in CCD readout. 

We used SExtractor to detect SSFs by looking for objects whose half-light radius is less than 1 pixel. We did not let 
SExtractor filter the images with matched-PSFs, and thus most of these SSFs would have occupied ~ 1 pixel if there 
had been no CTI-trails. We also measured the ellipticity of the SSFs by evaluating the unweighted moments using 
the unfiltered images. This ellipticity measurement scheme is different from Jee et al. (2009), where we measured 
the ellipticity after the image was convolved with a Gaussian function matching the instrument PSF. The qualitative 
behavior of the CTI does not change by this difference in measurement scheme. However, the current measurement 
more clearly reveals the "turnaround" feature that we discuss below. 

As a case study, we display the ellipticity bias of the SSFs in the galaxy cluster XMMXCS J2215-1738 in Figure [36} 
The left and right panels show the results for the SSFs, whose fluxes (counts) are high and low, respectively. Although 
here we use the SSFs detected only in the images of XMMXCS J2215-1738, they are still numerous and thus enable us 
to characterize the behavior of the CTI with high significance. First, it is clear that the average ellipticity of the SSFs 
in each flux range linearly increases (thus more negative) with the number of charge transfer. Second, for the bright 
SSFs, the CTI effect worsens for decreasing fluxes. In contrast, this trend is reversed for low-flux SSFs. This CTI 
turnaround was first reported in Jee et al. (2009) through both the SSF analysis and stellar photometry. Schrabback 
et al. (2010) also observed this turnaround in their independent CTI analysis. It is worth noting that the XMMXCS 
J2215-1738 cluster was observed in the year 2006 in short (^ 500 s) exposures, and the background level is very low 
(~ counts), which explains the severe CTI-trails. 

One way to compress the informati on in the two panels in Figure [36] is to plot the slopes of the data points as 



a function of flux as shown in Figure 31 The plot helps us to realize how sensitive the CTI slope is to the flux. 
Especially, it is remarkable that the CTI turnaround is very sharp (compare this with Figure 4 of Jee et al. [2009]). 

Before we proceed further, it is instructive to estimate what magnitude range of galaxies is mostly affected by this 
CTI pattern. As is mentioned in |2.2[ we select 24-28 mag galaxies as source objects. Considering that the average 
exposure time is ^-^500 s and the mean FWHM of the source galaxies is significantly larger than the SSFs (~ 8 pixels 
versus ~ 1.7 pixels), we estimate that this magnitude range corresponds to 2-100 counts in Figure 37 Therefore, the 
brightest and farthest (from the readout register) galaxies in our source sample suffer the CTI degradation most. 

B.2 CTI Correction to Galaxies 

We characterize the CTI with the ellipticity bias of the SSFs. Then, a nontrivial question is how to quantitatively 
relate this ellipticity bias to galaxy ellipticity change. Obviously, the ellipticity bias in SSF is a significant exaggeration 
of what occurs to galaxy shapes; a small CTI-trail induces much larger ellipticity to a 5 function-like feature than 
to, for example, a 0.5" galaxy. The exact translation requires our knowledge on the shape of the CTI-induced trail 
and the surface brightness profile of galaxies. Consequently, the recovery of the galaxy shape prior to the CTI effect 
involves the same delicacies in PSF-effect correction, which is a more familiar problem. 

Indeed, the CTI trailing can be treated as convolution with a one-dimensional kernel, which further smears the 
post-seeing objects, but only along the readout direction. 

o{x, y) = i{x, y) (g) p{x, y) (g) c{x, y) = i{x, y) ® r{x, y) (Bl) 

where i{x,y) and o{x,y) are the intrinsic and observed images, respectively. p{x,y) and c{x,y) are the convolution 
kernel representing the PSF and the CTI-trailing, respectively. Because the associative law holds for convolution, the 

When the slope is multipUed by 2048, the result gives the ■ n- i- i u- t iu i r ti, j i ■ j- 

^ 1' J 1 b maximum ellipticity bias tarthest irom the readout register 
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smearing of the object can be viewed as a single convolution by the kernel r{x,y), which includes both the PSF and 
CTI-trailing. One caveat here is that the CTI kernel c{x, y) in practice depends on flux, and, because a galaxy image 
consists of multiple pixels of varying flux, it is essential to approximate the effect with a single kernel, representing the 
collective effect of CTI on each galaxy. 

One can determine the shape of the kernel c(a;, y) by stacking many CTI-trails. We find that on average at short 
distances, CTI trails are well-approximated by an exponential tail cx e"^/"^, where the time constant t (expressed 
in units of pixel) determines how quickly the trail drops. Upon close examination, however, we observe that the 
measured CTI-trails differ in that 1) at large distances (> a few pixels) the measured CTI-trail slope is shallower, and 
2) occasionally a couple of small bumps appear, which suggests that different species of charge traps might be present. 
Nevertheless, we chose the simple functional form of a single exponential because object shapes are most sensitive to 
the information in the central few pixels. Based on the laboratory experiments of Dawson et al. (2008), Rhodes et 
al. (2010) also concluded that object shapes are most sensitive to the first few pixels and a single exponential term is 
adequate to describe the trails. 

Now the remaining question is how to determine the relation between the ellipticity bias 5e^ in the SSFs and the CTI 
time constant r. For this investigation, we rely on image simulations, where we create many artificial SSFs features 
and trail them with the ~ e'^l'^ kernel. Because we select only the ~ 1 SSFs to perform the study above, no 
exhaustive efforts are required to match the artificial SSFs to the real ones. Figure [38] shows the t vs. SSF ellipticity 
bias calibration that we derived from this image simulation. The thick solid line represents the relation when the 
size of the simulated SSFs matches the observed ones. The other lines illustrate how much the slope of the CTI trail 
depends on object sizes. The flowchart in Figure 39 summarizes the procedure to correct the PSF and CTI effects. 

Although the CTI correction described here is complicated and the result of time-consuming effort, we e mp hasize 
that the effect is on average small (also perhaps negligible in many cases) for cluster lensing analysis. Figure 40 shows 
that very few galaxies are subject to large elongation (> 0.05). A majority of galaxies, which fall to the CTI- mitigation 
regime are elongated by (5e+ << 0.01. 
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Fig. 30. — Point Spread Function correction. The left panels show the observed ellipticity pattern of the stars whereas in the middle 
panels we display the reconstructed ellipticity pattern at the location of the stars. The right panels compares the two ellipticity components 
e+ and ex before and after the correction. 
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Fig. 32. — Point Spread Function correction - continued from Figure [31] 
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Fig. 33. — Point Spread Function correction - continued from Figure [32] 
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Fig. 35. — Point Spread Function correction - continued from Figure [34] 
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Fig. 36.— CTI trail of SSFs in the XMIvIXCS J2215-1738 ACS images. The ellipticity bias (i.e., elongation along the parallel readout 
direction) due to CTI linearly increases (becomes more negative) with number of charge transfers. We divide the cases into two groups 
according to their fluxes in order to demonstrate the different flux-dependence. For the pixels with high counts (left), we observe that the 
CTI worsens for decreasing fluxes, as was reported in previous studies. However, the trend is reversed for the pixels with low counts (right). 
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Fig. 37. — CTI slope versus flux. Here we compress the information in the two panels in Figure [36] into one by plotting the slopes as a 
function of the fluxes. It appears that the flux-dependence can be classified into three regions. In the high-flux regime (> 500), the slope 
changes very slowly with flux. In the intermediate-flux regime (250 <counts< 500), there shows a precipitous decline of the CTI slope 
(more elongation) for decreasing flux. Finally, in the low-flux regime (< 250), a sharp turnaround is observed. We found that most of our 
background galaxies belong to the low end of this regime. Hence, the CTI effect is negligible for our source population. 




Fig. 38. — Time constant versus elongation. We model the CTI trail with oc e^^, where the time constant r determines how fast the 
trail truncates. The resulting ellipticity also depends on the size of the SSFs; for example, more ellipticity is induced for cosmic rays that 
occupy less number of pixels on average. The difference lines represent the results for different sizes. The thick solid line shows the case 
for the FWHM value of 1.65, which matches the mean size of the SSFs that we use for the study. 
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Fig. 39. — Galaxy shape measurement procedure with CTI and PSF corrections. 1. Surface brightness and coordinate are retrieved from 
the master source catalog. 2. Coordinate transformation is performed to locate the object in each visit's detector coordinate system. 3. 
Using the PSF model for each visit, we estimate the PSF at the object's location. 4. IBased on the y-coordinate (i.e., number of charge 
transfer pixels) and the surface brightness, we determine the ellipticity bias (e. g., ti sing the relation in Figure [371 1. Then, this ellipticity bias 
5e+ is converted to the CTI time constant t (e.g., using the relation in Figure [381 1. 5. The effective PSF (modified with CTI information) 
is rotated and scaled according to the exposure time prior to stacking for the creation of the final PSF. 
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Fig. 40. — Number of galaxies affected by CTI elongation. We examine how many background galaxies in the ACS field are affected by 
a given amount of CTI elongation (5e+. Very few galaxies are subject to large elongation (> 0.05). A majority of galaxies, which fall to the 
CTI-mitigation regime are elongated by (5e+ << 0.01, which can be safely ignored for cluster lensing analysis. 



